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Abstract 

We have constructed a Lagrangian model with a couphng of a and w mesons 
in the relativistic mean-field theory. Properties of finite nuclei and nuclear matter 
are explored with the new Lagrangian model SIG-OM. The study shows that an 
excellent description of binding energies and charge radii of nuclei over a large 
range of isospin is achieved with SIG-OM. With an incompressibility of nuclear 
matter K = 265 MeV, it is also able to describe the breathing-mode isoscalar 
giant monopole resonance energies appropriately. It is shown that the high-density 
behaviour of the equation of state of nuclear and neutron matter with the a-uj 
coupling is much softer than that of the non-linear scalar coupling model. 
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1 Introduction 



The relativistic mean-field (RMF) theory [T1[2l[3ll^ has been a successful ap- 
proach to describing properties of nuclei along the stability line as well as 
far away from it [5|l6|7j . Due to the Dirac-Lorentz structure of spin-orbit in- 
teraction, it has been shown to be advantageous over the conventional non- 
relativistic Skyrme theory in describing properties such as anomalous isotope 
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shifts in Pb nuclei [8]. It was shown [9] that an isospin dependence of the spin- 
orbit interaction or essentially a lack of it is responsible for the anomalous 
behaviour of the isotope shifts. The non-relativistic approaches based upon 
the Skyrme and Gogny forces have been unable to reproduce the kink [TOj . 
On the other hand, by significant alterations in the isospin dependence of the 
spin-orbit potential in the Skyrme theory pTj, it becomes possible to describe 
the above anomaly in isotope shifts. 

The prevalent model within the RMF theory is that of non-linear couplings 
of a meson [12|. Herein, the nonlinear scalar self-couplings of a meson play a 
pivotal role in describing finite nuclei. It was realized earlier on [12] that + 
(T^ terms have ingredients appropriate for a proper description of the surface 
properties of nuclei. By virtue of this, the scalar self-coupling terms have be- 
come indispensable for finite nuclei. One of the first nuclear forces constructed 
for calculation of ground-state properties of finite nuclei was NLl [2J. A larger 
value of the asymmetry energy of NLl, however, made its application to nuclei 
far away from the stability line difficult. Within this model, one of the first 
successful nuclear forces was NL-SH [6] which was obtained with a significant 
improvement in the asymmetry energy. The force NL-SH has an incompress- 
ibility of nuclear matter which is on the higher side of physical acceptable 
region. Sunsequently, force NL3 was obtained [7] with an improvement in the 
nuclear incompressibility of nuclear matter. This has in meanwhile become 
successful in several respects. 

A problem that afflicts the model with + cr^, in general, is that it gives 
an equation of state (EOS) of nuclear and neutron matter that is very stiff. 
Consequently, it is considered as incompatible with the observed spectrum of 
neutron star masses. A quartic term o;^ was added recently to the above model 
that was intended to improve the shell effects in nuclei [T^. It has subsequently 
been shown that the addition of the vector self-coupling of u meson has the 
advantage of softening the EOS of nuclear matter P^ffT3] . 

Variations in the relativistic approach based upon point- coupling model [T6|17fl8l 




or of density- dependent meson couplings [T9|20ll21ll22f23j have been made. 

With the inclusion of a larger number of parameters a good degree of success 
has been achieved [22]. There has also been significant discussion on the use 
of the effective field theoretical approach based upon expansion of the interac- 
tion Lagrangian in higher order terms in fields [2^25|26j . Addition of higher 
order terms and cross terms between the fields are shown to improve descrip- 
tion of properties of finite nuclei [2S]. Recently, a coupling between a and p 
meson has shown to be promising in reduction of neutron skin to reasonably 
accepted values as well as for obtaining a softer EOS of nuclear matter [27]. 
Thus, within the effective Lagrangian approach, the field of the RMF theory 
has become widely open with many possibilities of interaction terms. It is, 
however, not known as to what form of density dependence of nucleon-meson 
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interactions or of meson-meson interactions would be the most suitable one. In 
other words, which terms nature prefers in the expansion scheme is not known 
a priori. Emergence of nuclear data in extreme regions of the periodic table 
including those on finite nuclei, giant resonances and neutron star masses put 
nuclear interactions to test. 

With the initial success of the RMF theory in nuclear structure, there is thus 
a strong need to construct a theory that should be appropriate for broader 
aspects of finite nuclei as well as for nuclear matter. As improved and refined 
predictions of nuclear properties especially the masses (binding energies) of 
nuclei are required for regions often very far away from the stability line such 
as those for r-process nucleosynthesis, it is pertinent to devise new models and 
approaches to deliver the same. This may include other degrees of freedom or 
alternative forms of interactions. 

In this work, we have explored the RMF Lagrangian in the effective Lagrangian 
approach with a limited number of terms in the expansion with a view to see 
as to whether such an approach is feasible. We have restricted ourselves to a 
meson-meson interaction term between a and u mesons. Specifically, we have 
added a coupling of the form a^o;^, in addition to the + a'^ couplings which 
have so far been the most important terms in the RMF Lagrangian and are 
indispensable for finite nuclei and nuclear matter. 



2 Relativistic Mean-Field Theory: Formalism 

The basic (linear) RMF Lagrangian that describes nucleons as Dirac spinors 
interacting with the meson fields is given by [1] 



where is the bare nucleon mass and ip is its Dirac spinor. Nucleons interact 
with a, a;, and p mesons, where g^j, g^j, and gp are the respective coupling 
constants. The photonic field is represented by the electromagnetic vector A^. 

The nonlinear cr-meson-couplings added by Boguta and Bodmer ^12j are of 
the form 



Cq = ^{^^- g^i/j - gppf - -e(l - rg)^ - g^a - M^j ^ 
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UNL^lg2(T^+\93(T''. (2) 

The parameters g2 and gs are the nonhnear couphngs of cr-meson in the con- 
ventional + (j^ modeL The effective Lagrangian that is commonly used 
is 



Ceff = ^0 + Unl (3) 

Here, wc add an interaction term between the a and ou meson. Then, the 
effective Lagrangian becomes 



Ceff = jO,o + Unl + ^gatjcr'^i^ni^^- (4) 

The last term represents the coupling between a and uj meson that we in- 
troduce. The constant gcrw is the coupling constant of a and uj meson-meson 
interaction. The field tensors of the vector mesons and of the electromagnetic 
field take the following form: 



(5) 



The variational principle gives rise to the Dirac equation: 

{-ia.V + V{r)+f3m*}ipi^ e^^i, (6) 

where V{r) represents the vector potential: 

V{r) = g^cuoir) + g^npoir) + ei^Ao(r), (7) 

and S{r) is the scalar potential 

S{r) = g,a{v) (8) 

which defines the effective mass as given by 

m*(r) = m + ,S(r). (9) 

The corresponding Klein-Gordon equations can be written as 
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{-A + mf)a = -g^^Pij (10) 

(-A + m*2)^, = ^^^7,^/; (11) 

{-/\ + m^^p, = gp^l)^,TiP (12) 

-AA, = ]^eij{l + TM. (13) 

where, 

mf =ml + + g^a'^ - g^^uj^^ (14) 

*2 2 I 2 /I 



These equations imply an imphcit density dependence of a and uj meson 
masses. These density dependences are then responsible for the density de- 
pendence of nuclear interaction in a nucleus. This is in contrast to the non- 
relativistic Skyrme approach wherein the density functional is well-defined at 
the outset. An explicit density dependence of mesons has been introduced in 
refs. pU|22f23j . This approach requires several additional parameters besides 
the usual coupling constants in order to model a density dependence whose 
form is not known a priori. 

For the case of an even-even nucleus with a time-reversal symmetry, the spatial 
components of the vector fields u}, ps and A vanish. The Klein-Gordon equa- 
tions for the meson fields are then time-independent inhomogeneous equations 
with the nucleon densities as sources: 

(-A + ml)(j{r) = -gaps{v) - g2(y'^{v) - g^a^{Y) + g„^a{v)ujl{T) 
(-A + ml)ujQ{v) = g^p^{r) - g^^cx^{r)uJo{r) 

(16) 

(-A + m^)po(r) = gpPsir) 
-AAo(r) = ep,(r). 



For the mean-field, the nucleon spinors provide the corresponding source 
terms: 

A _ 

ps= E A 

i=l 
A 

pv= E i^t A 

N (17) 

P3 = E - E i^n 

'P=l n=l 

Z 

Pc= E V^p ^p, 
p=l 



where the sums are taken over the valence nucleons only. 
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The stationary state solutions ■0j are obtained from the coupled system of 
Dirac and Klein-Gordon equations self-consistently. The solution of the Dirac 
equation is obtained by using the method of oscillator expansion [5]. The 
ground-state of the nucleus is described by a Slater determinant |$ > of 
single-particle spinors -0, (i = 1,2,.. ..A). 

The centre-of-mass correction is added using the harmonic oscillator estimate: 
i?,,™, = '^-AlA-^'\ (18) 



3 The Lagrangian Parameter set 

The parameters of the new Lagrangian model are obtained by fitting bind- 
ing energies and charge radii of a set of nuclei within the RMF theory. The 
procedure has been described in ref. [6]. The nuclei included are ^^O, ^°Ca, 
48Ca, 9°Zr, ii^Sn, ^"^^^n, ^^ssn, and ^ospb. The Sn isotopes i^'^Sn, ^'^^Sii and 
^^^Sn are included with a view to take into account the broad range of isospin 
dependence of the nuclear interaction. For open-shell nuclei, the pairing is in- 
cluded within the BCS scheme, where the pairing gaps are obtained from the 
experimental masses of neighbouring nuclei. 

The binding energies and charge radii of nuclei have been used as constraints. 
In addition, we have tagged the spin-orbit splitting of — Pi/2 in ^^O. 
However, we have avoided putting explicit conditions or constraints on the 
nuclear matter properties. The u and p meson masses have been fixed at their 
empirical values. The procedure of obtaining the parameter sets of the new 
Lagrangian is the same as that applied in obtaining other forces such as NL-SH 
and NL3. It may recalled that the set NL3 was obtained by putting constraints 
explicitly on the nuclear matter properties as well. 

The parameters of the Lagrangian thus obtained (here named as SIG-OM) 
are listed in Table 1. These are compared with those of the forces NL-SH 
and NL3 with the nonlinear scalar self-couplings. The coupling constant 
with SIG-OM is positive in contrast with that of NL-SH and NL3, where it 
is predominantly negative. It has been argued [2] that a negative has a 
consequence in that the spectrum of the full theory is not bound from below 
and that renormalization of the scalar field is not possible. 

Solving the equations for nuclear matter at the saturation point, equilibrium 
properties of the nuclear matter are calculated. Nuclear properties arising from 
the parameters of SIG-OM are shown in Table 2. A comparison is made with 
the properties of NL-SH and NL3. The saturation density and binding energy 
per nucleon for SIG-OM are very close to those of NL-SH and NL3 and are in 
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Table 1 

The parameters of the Lagrangian SIG-OM with the couphng between a and u 
mesons. The Lagrangian sets with the non-hnear scalar self-coupling NL-SH, NL3 
are shown for comparison. 



Parameters 


SIG-OM 


NL-SH 


NL3 


M (MeV) 


939.0 


939.0 


939.0 


(MeV) 


505.9263 


526.0592 


508.1941 


ruu; (MeV) 


783.0 


783.0 


782.501 


nip (MeV) 


763.0 


763.0 


763.0 


9a 


10.0429 


10.4436 


10.2169 


9oJ 


12.7668 


12.9451 


12.8675 


9p 


4.4752 


4.3828 


4.4744 


92 (fm-1) 


-7.9223 


-6.9099 


-10.4307 


93 


12.4601 


-15.8337 


-28.8851 


9cru! 


35.6922 


0.0 


0.0 



Table 2 

The nuclear matter (NM) properties due to the force SIG-OM. Properties of the 
sets NL-SH and NL3 are also given for comparison. 



NM Properties 


SIG-OM 


NL-SH 


NL3 


PO (fm-=^) 


0.149 


0.146 


0.148 


ay (MeV) 


-16.30 


-16.33 


-16.24 


K (MeV) 


265.2 


354.9 


271.6 


m* 


0.622 


0.597 


0.595 


04 (MeV) 


37.0 


36.1 


37.4 



physically acceptable region. The incompressibility of nuclear matter K that 
arises from SIG-OM is 265.2 MeV. It is slightly smaller than that of NL3 
{K = 271.6 MeV). The effective mass at the saturation point with SIG-OM 
is m* — 0.62 . It is shghtly bigger than a value of ~0.60 for NL-SH and NL3. 
The asymmetry energy 04 of SIG-OM is 37.0 MeV. It is comparable to that 
of NL-SH and NL3 and is, however, still bigger than the received empirical 
value of -33 MeV. 
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4 Results and discussion 



4-1 Ground-state properties of spherical nuclei 



The binding energies and charge radii of nuclei calculated with SIG-OM in the 
RMF theory are presented in Tables 3. Results with NL-SH and NL3 are also 
given for comparison. Table 3 shows the binding energies achieved for the key 
nuclei included in the fit along with those of several other isotopes of Sn and 
Pb. The experimental binding energies [28] and charge radii where available 
are shown in the last column. The binding energies of the nuclei included in the 
fit are reproduced well by SIG-OM with a slight overbinding for ^^O. The set 
includes three doubly magic nuclei of ^°°Sn, ^^^Sn and ^°^Pb. It is important 
that any newly developed Lagrangian model should be able to describe the 
binding energies of closed-shell nuclei. Comparing the binding energy of ^^^Sn 
and ^''^Pb due to SIG-OM with the experimental data, one can see that there 
is a good agreement for ^'^^Sn. For ^°^Pb, SIG-OM overestimates the energy 
by ~ 1.5 MeV. This is, however, much improved as compared to that of NL- 
SH and NL3, where these sets overestimate the binding energy of ^osp]-, ]-,y 
~ 3 — 4 MeV. The binding energies of other isotopes of Sn and Pb nuclei 
obtained with SIG-OM show an excellent agreement with the experimental 
data as compared to NL-SH and NL3. 

The case of ^°°Sn needs a special mention. Almost all Lagrangian sets in use 
in the RMF theory overestimate the binding energy of ^''^Sn significantly. 
This includes both NL-SH and NL3 which overestimate the binding energy by 
~ 4 — 5 MeV. The prediction of SIG-OM for ^°°Sn is surprisingly close to the 
experimental value. This has, however, not affected (underpushed) the binding 
energies of ^^^Sn, ^^^Sn and ^^^Sn in the isotopic chain of Sn. Thus, SIG-OM 
provides the best available description of doubly magic nucleus ^°°Sn. 

The charge radii of nuclei calculated with SIG-OM are given in parentheses 
in Table 3. The SIG-OM values describes the experimental charge radii well. 
This is especially the case for Pb isotopes where SIG-OM values are closer 
to the experimental data. In comparison, NL3 overestimates the experimental 
charge radii of Pb isotopes. For the Sn isotopes, NL3 values are, however, closer 
to the experimental data. In comparison, charge radii for heavier isotopes 
with NL-SH are in close proximity to the experimental data. With the new 
Lagrangian model SIG-OM, there is a general agreement of its charge radii 
with the experimental data especially for heavier nuclei. 
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Table 3 

The binding energies in MeV and charge radii (parentheses) in fm of spherical 
nuclei obtained with SIG-OM. The values for NL3 and NL-SH are also shown for 
comparison. The empirical values (exp.) are shown in the last column 



Nucleus 


SIG-OM 


NL-SH 


NL3 


exp. 


160 


-129.2 (2.699) 


-128.4 (2.699) 


-128.8 (2.728) 


-127.6 (2.730) 


40Ca 


-343.3 (3.440) 


-340.1 (3.452) 


-342.0 (3.469) 


-342.1 (3.450) 


48Ca 


-414.6 (3.454) 


-415.0 (3.462) 


-415.1 (3.471) 


-416.0 (3.451) 


76Ni 


-631.8 (3.920) 


-634.4 (3.920) 


-634.1 (3.926) 


-633.1 


902r 


-783.1 (4.282) 


-782.9 (4.282) 


-782.6 (4.287) 


-783.9 (4.258) 


lOOg^ 


-826.0 (4.465) 


-830.6 (4.467) 


-829.2 (4.473) 


-824.5 




-988.8 (4.593) 


-987.9 (4.599) 


-987.7 (4.611) 


-988.7 (4.626) 


i24Sn 


-1049.9 (4.645) 


-1050.1 (4.651) 


-1050.2 (4.661) 


-1050.0 (4.673) 


i32Sn 


-1102.9 (4.698) 


-1105.9 (4.702) 


-1105.4 (4.708) 


-1102.9 


202pb 


-1591.8 (5.480) 


-1596.0 (5.480) 


-1592.6 (5.497) 


-1592.2 (5.473) 


208pb 


-1638.3 (5.506) 


-1640.4 (5.509) 


-1639.6 (5.523) 


-1636.7 (5.503) 


214pb 


-1663.8 (5.563) 


-1664.3 (5.562) 


-1661.6 (5.581) 


-1663.3 (5.558) 



4-2 The isotopic chain of Sn nuclei 



With the new Lagrangian model SIG-OM, we have performed study. 
The isotopic chain of Sn nuclei with the experimental masses available from the 
proton drip-line doubly-magic nucleus ^°°Sn to another doubly-magic nucleus 
^^^Sn offers unique data for a complete coverage of the shell from iV = 50 to 
= 82. We have performed the Relativistic Hartree-Bogoliubov calculations 
for Sn isotopes. The details of the method are provided in ref . [13] . The Gogny 
force DIS has been used in the pairing channel. It has been shown that the 
Gogny force is able to represent the pairing properties of Sn, Pb and other 
nuclei successfully PU] . 

The binding energies of Sn nuclei obtained with SIG-OM are shown in Fig. 1. 
A comparison is made with the predictions of NL3. Here, we have also included 
the results due to NL-SVl [13] • The force NL-SVl is based upon the model 
with the quartic coupling of u meson. The results show a significant difference 
between the predictions of SIG-OM and both the other Lagrangian models in 
the regions near the closed shells. Whereas NL-SVl shows a good agreement 
with the data in the region 106 < A < 132, NL3 shows a reasonably good 
agreement with the data in the region 106 < A < 126. Both NL3 and NL-SVl 
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Fig. 1. The binding energies of Sn isotopes calculated using the Relativistic 
Hartree-Bogoliubov method with the Lagrangian SIG-OM. The results obtained 
with NL3 and NL-SVl are also shown for comparison. 

show significant deviations from the experimental values below A < 106. For 
the double-magic nucleus ^°°Sn, the overbinding amounts to ~ 4 — 5 MeV. This 
problem of overbinding of magic nuclei and by implication the presence of a 
strong shell energy is not limited to theoretical models alone. The phenomenon 
of arches at the magic numbers persists both in microscopic theories and mass 
formulae alike. 

Viewing the binding energy curves in Fig. 1, one can see clearly that SIG-OM 
shows an excellent agreement with the experimental data all over the range of 
the shell. It shows a much smaller deviation from the data at A = 100. The 
well-known arch-like pattern that is exhibited by NL-SVl and NL3 is reduced 
significantly with SIG-OM. Thus, SIG-OM provides a very good description 
of the binding energies of the Sn isotopes over the whole region between the 
two magic numbers. 



4-3 Charge radii and isotopic shifts of Pb nuclei 



The salient feature of the Pb chain is the presence of a characteristic kink in 
charge radii and isotopic shifts at the magic number = 82. This behaviour 
of charge radii was considered for long to be anomalous [10] . Thus, the charge 
radii and isotopic shifts for Pb chain have become a standard test bench for any 
new model. In order to view as what response the new model SIG-OM gives 
for charge radii, we have performed RMF+BCS calculations for the isotopes 
of Pb. The charge radii of Pb isotopes calculated with SIG-OM are shown in 
Fig. 2. A comparison is made with the charge radii obtained with NL-SH and 
NL3. The experimental values [30] are shown for a comparison. A kink across 
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A = 208 (A^ = 82) arises in all the curves. The SIG-OM values describe the 
experimental data as well as NL-SH does. Whereas both SIG-OM and NL-SH 
show some minor differences with the data for isotopes lighter than ^°^Pb, 
for the heavier isotopes SIG-OM values show a better agreement with the 
data. The data point for ^''^Pb is significant as charge radius of this nucleus 
has been measured using various approaches and its value (rc = 5.503 fm) is 
known with a very high precision. The charge radius of ^ospb with SIG-OM is 
5.506 fm. It is very close to the experimental value. The value with NL-SH is 
5.509 fm which is slightly larger than the experimental data. In comparison, 
the charge radius of ^osp]^ with NL3 is 5.523 fm. The results of Fig. 2 show 
that NL3 overestimates the charge radii of the Pb isotopes systematically. 




Fig. 2. The charge radii calculated with SIG-OM, NL-SH and NL3. Experimental 
values are also shown for comparison 

The isotope shifts Ar^ = r'^{A) — r^(208) of Pb nuclei have been measured with 
a significant precision using laser-beam atomic spectroscopy [30]. Whilst ^osp^ 
acts as a reference point in the measurement of Ar^, the actual value of the 
charge radius of ^osp]^ does enter into the values of the isotopes shifts obtained 
from the measurements. We show in Fig. 2(b) the isotopes shifts calculated 
using the charge radii from SIG-OM, NL-SH and NL3. The experimental data 
[30] are shown directly. The theoretical values match the experimental data 
very well irrespective of the Lagrangian model used. Even the larger charge 
radii predicted by NL3 cancel out in the calculation of Ar^. Apparently, the 
isotope shifts of SIG-OM are very close to those due to NL3. We can con- 
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elude that the RMF theory describes the experimental isotope shifts very 
well irrespective of the nature of the Lagrangian model employed with barely 
perceptible differences amongst various forces. 



4.4 Nuclei away from the stability line 

Mass measurements with high precision on nuclei far away from the stabil- 
ity are being performed at various facilities and new experimental data are 
emerging currently. Such data provide a crucial test bench for various nuclear 
models in regions away from the stability line. In order to test the predictive 
power of the Lagrangian SIG-OM, we have performed deformed RMF calcula- 
tions for some nuclei on which experimental data has been obtained recently. 
The total binding energy of some isotopes of Si, Sr, and Mo and a few others 
obtained with SIG-OM is shown in Table 5. A comparison is made with NL3. 
The experimental binding energies obtained from high-precision experimental 
masses on Si [21], Sr [22]) and Mo [33]) are used for comparison. 



Table 4 

The binding energy E (in MeV) and quadrupole deformation P2 of nuclei away from 
the stability line calculated with SIG-OM. The deviation of the binding energy from 
the experimental value, AE = E - Eg^p, is also shown. A comparison is made with 
the predictions of NL3. 



Nucleus 




SIG-OM 






NL3 






E 


P2 


AE 


E 




AE 


36Si 


-292.4 


0.02 


-0.4 


-293.5 


0.0 


-1.5 


38Si 


-300.0 


0.28 


-0.2 


-301.3 


0.27 


-1.5 


40si 


-306.6 


0.36 


-0.1 


-307.9 


0.34 


-1.4 


42Si 


-313.3 


-0.35 


-0.4 


-315.2 


-0.33 


-2.3 


soSr 


-684.8 


0.0 


+1.5 


-682.8 


0.0 


+3.5 


s^Sr 


-748.9 


0.0 


0.0 


-748.2 


0.0 


+0.7 


88Sr 


-768.1 


0.0 


+0.4 


-768.0 


0.0 


+0.5 


108M0 


-909.2 


-0.22 


+0.4 


-908.3 


-0.23 


+1.3 


"OMo 


-919.6 


-0.23 


-0.1 


-919.0 


-0.23 


+0.5 


i20Xe 


-1008.2 


0.26 


+0.3 


-1007.7 


0.29 


+0.8 


174Yb 


-1407.2 


0.31 


-0.6 


-1407.2 


0.32 


-0.6 
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It is noted that SIG-OM describes the experimental binding energies very 
well with only a few exceptions. The agreement is especially good for newly 
measured Si isotopes. Not surprisingly the force NL3 also shows reasonably 
good agreement for several data points. It, however, shows somewhat larger 
deviations from the data notably for the Si isotopes. Details of a comprehensive 
study to examine the shell structure and shell effects due to the Lagrangian 
model SIG-OM encompassing a larger number of data in the regions far away 
from the stability line will be discussed elsewhere |34] . 

4.5 Breathing mode GMR energies and GCM calculations 

The incompressibility of nuclear matter K is an important point on the equa- 
tion of state (EOS) of nuclear matter. By definition, it is the curvature of the 
EOS at the saturation point. The value of K is important as it affects the 
breathing mode giant monopole resonance (GMR) energies sensitively. As the 
force NL-SH has a rather large value, it is found to be not suitable for the 
GMR energies of nuclei. It may be recalled that within the + model, 
the force NL3 was obtained with a value of incompressibility K = 271 MeV 
which is in the range of physically acceptable values. In order to get a compar- 
ative picture of the response of various Lagrangians on the breathing-mode 
giant monopole resonance (GMR) energies, we have performed constrained 
generator coordinate method (GCM) calculations for the GMR mode [35]. 
The nuclei included are ^°Zr, ^^°Sn and ^°^Pb. Experimental GMR energies on 
these nuclei are known with a reasonable precision [36|l37] . The results of our 
calculations are shown in Table 6 and are compared with the experimental 
data. Clearly, the force NL-SH [K = 355 MeV) provides larger GMR energies 
due to its high incompressibility. NL3 {K = 271 MeV), on the other hand, un- 
derestimates the GMR data by ~ 1 MeV. Results of our calculations obtained 
with NL3 agree well with those provided in ref. [7]. 

Table 5 

The breathing mode GMR energies in nuclei obtained with the the constrained 
GCM calculations. The experimental data are from refs. [36|37] 



Nucleus 


SIG-OM 


NL3 


NL-SH 


exp. 


90Zr 


18.2 


16.9 


19.5 


17.81 ±0.30 


i20Sn 


16.2 


15.0 


16.7 


15.52 ±0.15 


208pb 


14.1 


13.0 


15.0 


13.96 ± 0.28 



It is interesting to note that SIG-OM with its slightly lower value of i^' = 265 



13 



MeV than that of NL3 provides GMR energies which are closer to the ex- 
perimental data, especially those of ^^Zr and ^°^Pb considered frequently in 
analyses. In comparison, NL3 with a higher value oi K = 271 MeV under- 
estimates the data by ~1 MeV. Generally, a force with a lower value of K 
would yield smaller values of GMR energies. This difference in the predictions 
of SIG-OM and NL3 may be due to a difference in some finite size effect of 
a nucleus. The other two factors which play a role for GMR energies are the 
Coulomb term and the asymmetry component. At present, it is difficult to 
estimate as to which factor brings about the apparent paradoxical behaviour 
of SIG-OM and NL3. We believe that it may be due to a different surface 
component (arising out of a different density dependence of the interactions) 
in SIG-OM vis-a-vis NL3. It will, therefore, be important to carry out calcula- 
tions of surface incompressibility using the semi-infinite nuclear matter. These 
calculations are in progress. Moreover, it is also to be seen as to what response 
the two Lagrangian models provide in relativistic radom-phase approximation 
calculations of the breathing-mode GMR energies. 



4-6 Nuclear matter properties 
4-6.1 Effective mass 

The density dependence of the effective mass m* is to an extent a reflection 
of the density dependence of interactions in a nucleus, though it corresponds 
to that of the a field to be exact. It would be instructive to see as to how 
m* changes with density in different Lagrangian models. We show in Fig. 3 
the density dependence of m* for different Lagrangian models. Fig. 3(a) shows 
the effective mass m* for the density region of a nucleus up until a 4-5 times 
the saturation density. Both NL-SH and NL3 bunch out together and show a 
small difference in the high density region. NL-SVl sphts out from NL-SH and 
NL3 at higher densities. In the high-density region, three forces NL-SH and 
NL3 with the scalar self-couplings and NL-SVl with the scalar self-coupling 
along with the vector self-coupling exhibit a similar pattern of variation with 
density as seen in Fig. 3(a). The model SIG-OM with a — lo couphng, on the 
other hand, shows a density dependence that is different from that of NL-SH, 
NL3 and NL-SVl and it is especially so in the high density region. 

The region of the density relevant to finite nuclei is shown in Fig. 3(b). The 
saturation point is indicated by a square in each case. Both SIG-OM and NL- 
SVl have a slightly larger m* than that of NL-SH and NL3. It may be pointed 
out though that in each case the spin orbit splitting of ^^O is reproduced well 
and yet overall fits to finite nuclei lead to slightly different effective mass for 
each model. This is no coincidence that within the + o"*^ Lagrangian model, 
the effective masses for NL-SH and NL3 are very close to each other. The 



14 




Fig. 3. The effective mass m* calculated with SIG-OM. It is compared with that 
for the 0"^ + cr^ Lagrangians NL-SH and NL3. A comparison is also made with that 
of (j^ Lagrangian NL-SVl. (a) The density dependence of m* for the whole density 
density including the high-density region, (b) m* in the region of density relevant 
to finite nuclei. The saturation point for each case is shown by squares. 

variation in m* for these two forces in this region is also very similar. On the 
other hand, in the region of interest for finite nuclei, it is seen clearly that 
SIG-OM shows a density dependence that deviates strongly from that of NL- 
SH and NL3. This may possibly lead to a different features of surface with 
SIG-OM. Detailed surface properties using various Lagrangian models are in 
progress. It is a matter for further investigation as to whether this difference in 
density dependence leads to different predictions of breathing-mode energies 
with SIG-OM as discussed above. 



4-6.2 The EOS of nuclear and neutron matter 

The EOS of nuclear and neutron matter is important for structure and proper- 
ties of neutron stars. There has recently been a significant discussion on EOS 
of nuclear matter |38] especially in view of recent observations of neutron stars 
which have been found to exist with masses in the vicinity of two solar masses 
or even more [SH]. In accordance with the known spectrum of neutron star 
masses which lie below ~ 1.8 solar mass (Mq), the conventional wisdom has 
been that a softer EOS is required for describing the observed masses. How- 
ever, recent discoveries of masses close to 2Mq or more than 2Mq seem to 
put new demands on the high-density EOS of nuclear matter. In ref. , the 



15 



neutron star EXO 0748-676 with mass M > 2.1O±O.28M0 has been reported. 
Consequently, very soft EOS for nuclear and neutron matter are ruled out. 
Then, one requires an EOS that should encompass even the heavier neutron 
star masses in the region of 2.0 — 2.3Mq. 

Solving the equations for nuclear matter self-consistently, we have calculated 
the EOS of nuclear and neutron matter with the new Lagrangian model. The 
EOS of nuclear matter with SIG-OM shows a significant softening at higher 
densities as compared to that of NL-SH and NL3 (see Fig. 4). The coupling 
of a and uj meson does have a softening effect at higher densities. In contrast, 
EOS with the + model is known to be stiff due to preponderance of a; 
term at higher densities. 



Equation of State of Nuclear Matter 



Neutron Matter EOS 




P (fm-1 



Fig. 4. The EOS of nuclear matter and neutron matter (inset) obtained with 
SIG-OM, NL-SH and NL3. A comparison is made with the EOS due to Friedman 
and Pandharipande (FP) [iO] . 



The corresponding EOS for neutron matter is shown in the inset of Fig. 4. 
It can be seen that neutron EOS due to SIG-OM is also softer than that for 
NL-SH and NL3. For a comparison, we show the EOS due to Friedman and 
Padharipande (FP) [40]. The FP EOS for nuclear as well as neutron matter is 
very soft and was intended to describe lower mass neutron stars with masses 
M < 1.6M0. As one can see from Fig. 4, the softening achieved by the a — u 
coupling is not strong enough to match that of FP which is a very soft EOS 
leading to smaller values of maximum mass of neutron stars [lO] . In view of 
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new observations such as those of ref. [39], a softer EOS in the high density 
region may not be required. It has to be seen from neutron-star structure 
calculations whether the EOS due to SIG-OM which is neither as stiff as NL3 
nor as soft as FP would be suitable for describing the newly observed neutron 
star masses. 



5 Conclusions and outlook 

The Lagrangian SIG-OM with the coupling of a and uj mesons in the RMF 
theory has been developed. We have explored the feasibility of adding the a-uj 
coupling in the RMF Lagrangian. Properties of finite nuclei and nuclear matter 
have been explored with the model SIG-OM. It is shown that SIG-OM provides 
a very good description of the ground-state properties such as binding energies 
and charge radii of nuclei over a large range of isospin. Taking the case of the 
isotopic chain of Sn nuclei, it is shown that it describes the binding energies 
of nuclei between the two magic numbers = 50 and = 82 very well. It 
will be interesting to see in further investigations as to what implications this 
improvement in shell structural aspects may have in the other regions and 
most importantly in the regions close to the r-process path. 

The breathing-mode giant monopole resonance energies for a few nuclei have 
been calculated using the constrained GCM calculations. It is shown that the 
model SIG-OM with the incompressibility of nuclear matter K = 265 MeV 
describes the breathing-mode GMR energies well. A comparative analysis of 
breathing-mode GMR energies with the GCM approach in the RMF theory 
suggests that some finite-size effects with SIG-OM may be different than with 
Lagrangians with nonlinear scalar self-couplings. It is important to carry out 
calculations with relativistic RPA approach to see as to whether such a feature 
would persist. This is a matter of further investigations and requires a detailed 
study. 

The equation of state of nuclear and neutron matter have been obtained with 
SIG-OM. A comparison is made with the EOS obtained with the interactions 
within the nonlinear scalar self-coupling. As it is well-known, the latter model 
gives very stiff EOS. The equation of state of nuclear and neutron matter 
obtained with the new Lagrangian SIG-OM is significantly softer than that 
with non-linear scalar self-couplings. 

We thank Prof. Lev Savushkin for fruitful discussions. 
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